A new measurement of the 2: lifetime 



The FOCUS Collaboration* 

J. M. Link a M. Reyes a P. M. Yager a J. C. Anjos b I. Bediaga b 
C. G6bel b J. Magnin b A. Massafferri b J. M. de Miranda b 

I. M. Pepe b A. C. dos Reis b S. Carrillo c E. Casimiro c 
E. Cuautle c A. Sanchez-Hernandez C. Uribe F. Vazquez c 

L. Agostino d L. Cinquini d J. P. Cumalat d B. O'Reilly d 
J. E. Ramirez d I. Segoni d J. N. Butler e H. W. K. Cheung 
G. Chiodini 6 I. Gaines 6 P. H. Garbincius L. A. Garren e 
E. Gottschalk P. H. Kasper A. E. Kreymer R. Kutschke 
L. Benussi f S. Bianco f F. L. Fabbri f A. Zallo f C. Cawlfield g 
D. Y. Kims K. S. Park§ A. Rahimis J. Wiss§ R. Gardner 11 
A. Kryemadhi h K. H. Chang 1 Y. S. Chung 1 J. S. Kang 1 
B. R. Ko 1 J. W. Kwak 1 K. B. Lee 1 K. Cho j H. Park^ 
G. Alimonti k S. Barberis k A. Cerutti k M. Boschini k 
P. D'Angelo k M. DiCorato k P. Dini k L. Edera k S. Erba k 
V ■ M. Giammarchi k P. Inzani k F. Leveraro k S. Malvezzi k 

<D . 
43 ■ 



X 



D. Menasce k M. Mezzadri k L. Moroni k D. Pedrini k 
C. Pontoglio k F. Prelz k M. Rovere k S. Sala k 
T. F. Davenport lll e V. Arena 111 G. Boca m G. Bonomi m 
G. Gianini m G. Liguori 111 M. M. Merlo 111 D. Pantea m 
S. P. Ratti m C. Riccardi 111 P. Vitulo m H. Hernandez 11 
A. M. Lopez 11 H. Mendez 11 L. Mendez 11 E. MontieP D. Olaya 11 

A. Paris 11 J. Quinones 11 C. Rivera 11 W. Xiong n Y. Zhang 11 
J. R. Wilson T. Handler p R. Mitchell p D. Engh q M. Hosack q 
W. E. Johns q M. Nehring q P. D. Sheldon q K. Stenson q 
E. W. Vaandering q M. Webster q M. SheafF 

a University of California, Davis, CA 95616 
h Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brasil 
C CINVESTAV, 07000 Mexico City, DF, Mexico 
d University of Colorado, Boulder, CO 80309 



see http://www-focus.fnal.gov/authors.html for additional author information. 



Preprint submitted to Elsevier Science 



7 February 2008 



e Fermi National Accelerator Laboratory, Batavia, IL 60510 
i Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy 1-00044 
g University of Illinois, Urbana- Champaign, IL 61801 
^Indiana University, Bloomington, IN 47405 
1 Korea University, Seoul, Korea 136-101 
i Kyungpook National University, Taegu, Korea 702-701 
^INFN and University of Milano, Milano, Italy 
£ University of North Carolina, Asheville, NC 28804 
111 Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy 
11 University of Puerto Rico, Mayaguez, PR 00681 
° University of South Carolina, Columbia, SC 29208 
p University of Tennessee, Knoxville, TN 37996 
q Vanderbilt University, Nashville, TN 37235 
1 University of Wisconsin, Madison, WI 53706 



Abstract 

Using data collected by the Fermilab experiment FOCUS, we measure the 
lifetime of the charmed baryon S° using the decay channels S° — > S~7r + 
and S° — > Q~K + . From a combined sample of 110 ± 17 events we find 
r(S°) = 118 1;^ =t 5 fs, where the first and second errors are statistical 
and systematic, respectively. 



1 Introduction 



While the lifetime hierarchy of the weakly decaying charm mesons is well 
established experimentally [1] and in fairly good agreement with theoretical 
calculations [2] , the pattern of lifetimes in the charm baryon sector still needs 
to be established both experimentally and theoretically [3]. While the A+ 
lifetime is known to an accuracy of about 2% [4-6], and recently the 
lifetime has been determined to an accuracy of 5% [7], the S° and lifetimes 
remain known to only 20% and 30%, respectively. 

The lifetime of the S° has been measured by two experiments. CERN experi- 
ment NA32 found a lifetime of 82 j^ fs using four events reconstructed in the 
decay channel S° — > pK~K*° [8]. Fermilab experiment E687 [9] measured a 
S° lifetime [10] of 101 jiff (stat) ±5 (sys) fs using 42 events from the decay 
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FOCUS (Fermilab E831) is a general purpose experiment investigating charm 
physics. The charmed particles are produced by the interaction of a photon 
beam on a segmented BeO target. The average energy of the beam for the 
data collected for this measurement is about 180 GeV. The track reconstruc- 
tion is accomplished by two silicon vertex detectors [11,12] (TS and SSD) 
in the target region and by five multi-wire proportional chambers (MWPC) 
downstream of the target region. The charged momenta are determined by a 
measurement of the bending angles in two magnets of opposite polarity. Due 
to the excellent separation between production and decay vertices provided 
by the silicon detectors, we achieve an average proper time resolution of about 
50 fs for the decay channels used in this analysis. Charged particle identifica- 
tion of hadrons is performed with three multi-cell threshold Cerenkov counters 
[13]. 



2 Reconstruction of hyperons S and f2 

We reconstruct the decay modes S° — >■ and S° — > Q~K + . A detailed 

description of the reconstruction method and the mass spectra of the hyperons 
2~ and is in Reference [14]. We reconstruct the S~ decay channel A°7r~ and 
the Q~ decay channel A°K~ . While these decay topologies are very similar, we 
can easily distinguish between them using Cerenkov identification of the final 
state particles and by requirements of the reconstructed invariant masses. We 
reconstruct the entire decay chain of these two modes. The A is identified by 
its decay into a proton and an oppositely charged pion. The tracks of these 
charged daughters are used to form the decay vertex of the A and to determine 
its flight direction and momentum. This direction is used together with the 
momentum of either the 7T~ or the K~, to determine the decay vertex and 
momentum of the hyperon. The hyperon vertex must lie upstream of the A 
vertex. Finally, the position and momentum vector of the hyperon is matched 
to a SSD track. 

For both modes we compute the invariant mass of the two body combination 
and select the events which are in the correct hyperon mass region. For the 
S" we require 1.312 GeV/c 2 < M(A°tt-) < 1.330 GeV/c 2 and for the Q~ we 
require 1.665 GeV/c 2 < M(A°K~) < 1.680 GeV/c 2 . 



3 Reconstruction of S° candidates 

The S° candidates are reconstructed using a candidate driven vertexing algo- 
rithm [9]. Briefly, combinations of tracks are used to form a secondary vertex 
which in turn must point within errors to a primary production vertex. Pairs 
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of H~ candidates and oppositely charged pions (kaons) are intersected to 
form a secondary vertex and a S° mass. Several criteria are used to reject the 
background contributions. Any track consistent with an e + e~ pair hypothesis 
is rejected. A confidence level (CLS) for the secondary vertex hypothesis from 
the intersection of the two tracks is calculated, and a minimum cut is applied, 
the magnitude of which depends on the decay mode. The S° decay products 
are used to construct the S° momentum vector, which is combined with at 
least two additional silicon tracks to form the primary vertex, which must lie 
within the target region. A confidence level (CLP) for the primary vertex is 
calculated, and the combination is rejected if CLP < 1%. 

We require a minimum value for the significance of the separation between 
the secondary and primary vertex, given by i/cre, where £ is the distance be- 
tween the two vertices and ae is the uncertainty on £. The Cerenkov particle 
identification uses requirements on the variables W e , W n , W K , and W p , which 
represent the negative log likelihood for the hypothesis that a particular track 
is an electron, a pion, a kaon, or a proton, respectively. The difference be- 
tween two of these variables represents the ratio between the two probabilities. 
The pion consistency of a track is defined by a requirement on the variable 
picon = W m i n — W n , where W m i n is the minimum likelihood of the four corre- 
sponding hypotheses. The specific requirements for the two decay modes are 
now discussed. 

For the mode, we apply the detachment cut i/ag > 3 and we require 

CLS > 3%. During the data collection the tracking system was improved with 
the introduction of a silicon vertex detector in the target region; we require 
o~ t < 150 fs and o~ t < 80 fs where o~ t is the error on the decay time of the 
S°, respectively, for data collected before and after the introduction of this 
detector. The pion (decay product of the S°) is required to have picon > —10 
and a momentum greater than 7 GeV/c. 

For the Q~K + mode, the detachment cut is tjo~i > 0.5. For the secondary 
vertex we require CLS > 1%. To identify kaons, we require W p — Wk > —0.1 
and W n — Wk > 1, comparing the kaon hypothesis to the proton and the 
pion hypothesis. Since S~ particles are produced more copiously than Q~ 
particles, we compute the invariant mass for the fl~ when it is calculated as 
a A°7r~ combination. To eliminate S~ contamination under the Q~ mass, we 
require M(A°tt-) > 1.375 GeV/c 2 . 

Figure 1 shows the invariant mass plots for the two modes. The fit is to a 
Gaussian for the signal and a straight line for the background. We do not use 
the fit information to perform the lifetime measurement. 
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4 The lifetime measurement 



We use a binned maximum likelihood fit method to measure the lifetime. The 
fitted histogram is the reduced proper time distribution: t' = (£ — Nae)/((3 / yc), 
where N is the vertex detachment cut, and ^7 = p(S°)/M(S°). The fit is 
performed on the data in the signal region of the invariant mass distribution: 
the mass window within ±2 a of the S° nominal mass, where a is the width 
of the Gaussian fit of the invariant mass distribution. The t' distribution of 
the signal region differs from the expected behavior e~ 1 ' l T due to the presence 
of background events. We assume that the background events in the signal 
region have the same time distribution as the events of the sidebands (located 
between 5 and 9 a away from the S° mass). If S is the total number of signal 
events in the signal region, and B is the total number of background events 
in the signal region, the expected number of events rii in the i th bin of the t' 
distribution is: 

n , =s mx, + B.. b - 



where bi is the number of events in the i th bin of the sideband t' distribu- 
tion. The Monte Carlo correction function f(t') takes into account detector 
acceptance, the efficiency of the selection cuts, and absorption of the daugh- 
ter particles as a function of the reduced proper time. Figure 2 shows the 
correction functions for the two modes; due to the lower detachment cut we 
observe a bigger correction for the Q~K + mode. Figure 3 shows the reduced 
proper time distribution for each mode after the Monte Carlo correction. For 
each mode a likelihood is constructed as the product between the Poisson 
probability of observing Sj events in the signal region, when rii are expected, 
and the Poisson probability of observing J2i h background events when B are 
expected. A factor of 2 is included because the sidebands are twice as wide as 
the signal region. The likelihood for each category is thus: 

L - {Ui ^r ) x ( <zM 1 

To measure r, we maximize the combined likelihood, which is the product of 
the likelihoods for the two decay modes: 

L~o = Lsn x Lq K . 

There are three fit parameters: the lifetime r, and for each mode the back- 
ground events B. Our fit result is r = 118li 2 When fitted separately, we 
measure for the two decay modes r(S° — > 2~7r + ) = 1171 i 3 fs and r(S° — > 
tt~K+) = 1221^ fs. 
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Contribution 


Systematic (fs) 


Production 


4- Q 
± o 


Fit 


± 4 


Method 


± 2 


Total 


± 5 



Table 1 ~ ~ ~ " 

Contributions to the systematic uncertainty. 

5 Systematic uncertainty determination 

The systematic uncertainty was evaluated from a detailed study of different 
sources: the Monte Carlo production and detector simulation, the fitting pro- 
cedure, and the possible intrinsic bias of the method. 

The Monte Carlo correction function may be incorrect if the Monte Carlo 
poorly simulates the detector or the production characteristics of the events. 
For this reason we extract a systematic uncertainty due to the Monte Carlo 
simulation after a careful study of how well the simulation matches the data. 
We analyze effects from the momentum and transverse momentum of S[j, the 
flavor (particle and antiparticle) of the the decay mode of the S°, the mul- 
tiplicity and position of the primary vertex, the error on the decay time and 
length, and the target region silicon strip detector (available for 2/3 of the 
FOCUS data set). To obtain a systematic error, we split the data into statis- 
tically independent samples, on the basis of the production variables and of 
the mode. Then we compute the x 2 /d.o.f. for the hypothesis of consistency of 
the measurements. If we find x 2 /d.o.f. > 1, we rescale the errors in order to 
have x 2 /d.o.f. = 1, and extract the systematic error subtracting in quadrature 
the statistical error from the scaled error of the weighted average of the inde- 
pendent measurements. This method is based on the ^-factor method of the 
Particle Data Group [1] . The only significant effects found are due to primary 
vertex multiplicity and S° decay mode from which a systematic uncertainty 
of cr(production) = ± 3 fs is obtained. 

The measurement we report is from a particular choice for the fit parameters 
and fitting technique (namely, the binned maximum likelihood technique). 
Varying these fitting conditions provides measurements which are all a priori 
likely. We calculate a systematic uncertainty by performing a set of lifetime 
measurements with different choices for the fitting parameters and the fitting 
technique. To study systematic effects related to our choice of fitting parame- 
ters we varied the location and width of the sidebands, the proper decay time 
fitting region and bin size. We also varied the correction function f(t'), by 
using a linear fit instead of the bin values. 



6 



Since the proper time resolution (about 40 fs for the 2~7r + mode and 80 fs 
for the Q~K + mode) is close to the S° lifetime, we decided to perform the 
measurement also with a convolved binned likelihood method [15] to check the 
fitting technique. In this method, the exponential decay is convolved with the 
smearing due to the time resolution. We used the Monte Carlo convolution cor- 
rection function F(t^ en , t\ ec ), which represents the probability of reconstructing 
an event in the decay time bin % when the true decay time bin is j. The like- 
lihood is constructed in the same way as explained before, but the number of 
expected events in the i th bin becomes: 



n, 



~ Ei Ej f(tf n , tf c )e-^/ T Y,ih 



Instead of the reduced proper decay time, the proper decay time is fit. Fig. 4 
shows the proper time distributions. The lifetime result of the convolved 
method is then considered as a further fit variant. The systematic uncertainty 
due to the fit variants is given by the variance of the set of measurements. We 
find a(fit) = ± 4 fs. 

The Monte Carlo simulation uses an input lifetime given by the current world 
average reported from the PDG [1] of t(S°) = 98 fs. We investigate the range 
of good performance of the fitting method when the input lifetime for the 
correction function is 100 fs with a mini Monte Carlo study. This consists of 
simulating a large number of samples with the same statistics as our data (both 
for signal and background events). We obtain events from exponential decaying 
populations with a given lifetime for the signal and for the background. We 
obtain the background lifetime by fitting the data sideband distributions to 
an exponential. For the signal events we study a wide range of lifetime values 
(from 10 fs to 170 fs). Care is taken to account for the time resolution. The 
method works with an accuracy better than 10% for lifetimes greater than 
55 fs. We study the accuracy in the range of the S° measured lifetime, and 
measure a systematic uncertainty a (method) of ± 2 fs. The mini Monte Carlo 
study also validates our estimate of the fit statistical error. 

The total systematic uncertainty is obtained by adding in quadrature the con- 
tributions from the three independent sources. We thus find a(sys) = ± 5 fs. 
See Table 1 for a summary of the contributions and the total systematic un- 
certainty. 



6 Conclusions 



We have measured the lifetime of the charmed baryon S°, fully reconstruct- 
ing the decays S° — > and S° — > £l~K + . Figure 5a shows the invariant 
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mass plot for the combined sample. From the 110 ± 17 reconstructed events 
we measure t(S°) = H8I12 (stat) ± 5 (sys) fs. Figure 5b shows the Monte 
Carlo corrected sideband subtracted proper time distribution for the com- 
bined sample for the signal region and sideband region. This measurement 
greatly improves upon the accuracy of previous measurements, reducing the 
percentage error on the S° lifetime from 20% to 10%. 
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Fig. 1. Invariant mass distributions for each mode. The fit uses a Gaussian for the 
signal and a straight line for the background. The signal region is defined by the 
central shaded region. The sideband region, defined by the two symmetrical shaded 
regions, contains the events used to reproduce the decay time of the background 
events in the signal region. 
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Fig. 2. Monte Carlo correction function f(t') for each decay mode. 




Fig. 3. Reduced proper time distributions, corrected by the Monte Carlo correction 
function, for each decay mode in the signal region and the sideband region (shaded). 
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Fig. 4. Proper time distribution for the two decay modes in the signal region and 
sideband region (shaded). 
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Fig. 5. Mass and lifetime distributions for the combined sample, a) Invariant mass 
distribution; fit with a Gaussian for the signal and a first order polynomial for the 
background, b) Monte Carlo corrected and sideband subtracted reduced proper time 
distribution for the signal region (points). The histograms with different shades are 
the time distributions for the sidebands for the two decay modes (lighter shade: 
H~7r + , darker shade: Q~K + ). The dashed line shows the result of the fit. 
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